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Abstract 

The fluctuating field of a jet excited by transient mass injection is 
simulated numerically. The model is developed by expanding the state vector 
as a mean state plus a fluctuating state. Nonlinear terms are not neglected 
and the effect of nonlinearity is studied. The results show a significant 
spectral broadening in the flow field due to the nonlinearity. In addition, 
large scale structures are broken down into smaller scales. 
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Introduction 


A numerical method for simulating the axisymmetrlc fluctuating field of a 
forced circular jet has been developed. The method is based on a solution of 
the Euler equations for homentropic flow in cylindrical coordinates. The jet 
is forced by transient point mass injection. _T^_squrce strength is scaled by 
a parameter e and the state vector W = (p,p,u,v) is expanded as a_ mean 
state plus a fluctuating state of order e. Here p is the pressure, p the 
density and u and v are the axial and radial velocity components respec- 
tively. 

The solution is expanded in the form p = p. + ep, etc. where the mean 
state (with subscript "0") is assumed to satisfy the unforced Euler 
equations. The variables (p,p,u,v) will represent the fluctuating field in 
response to the jet forcing. New variables e = (Pj.+ep)u and 
f = (Pg+ep)v are introduced for simplicity. With cylindrical spatial co- 
ordinates z (axial) and r (radial) the Euler equations for the fluctuating field 
become 

Pt + (pUg+Pgu)^ + (pVg+Pgv)^ + ^ 

®t te(Uo+eu))^ + (e(VQ+ev))^ + p^ + fU^^ - eV^^ - = emu (1) 

ft + (fCUg+Eu))^ + (f(VQ+ev))^ + Pr " ^^Oz " 

In the development of system (1) it is assumed that the derivatives of 
the mean pressure and density in space can be neglected. It is an exact con- 
sequence of the full Euler equations with a source of mass Injection, and of 
the expansion into mean and fluctuating states. The source is em (units 

density/time) . The mean velocities Uq and Vq are taken from measurements 
of a spreading jet. The system is solved in a cylindrical rectangle including 
a semi-infinite pipe from which the jet exits. The computational domain is 
Illustrated in Fig. 1. 

The use of the quasi-momentum variables e and f gives a system from 
which the fluctuating field can be computed directly rather than as a small 
part of the total field* In addition the linear limit can be recovered by 
simply setting e = 0 in (!)• The nonlinear terms are explicitly 
exhibited* The fluctuating pressure p is obtained from the density p by 
the homeintropic relation 

p = = Pq(1 + £ -^)^, (2) 

where y = in air. 

The fundamental assumption entering into the derivation of (1) is that 
the mean state is a solution to the unforced Euler equations. This is not 

true for a state determined from experimental measurements. Numerical experi- 
ments have verified however, that the qualitative features of the fluctuating 
solution are insensitive to small changes in the mean state and thus we 
believe that the solution to (1) qualitatively represents the fluctuating 
field in response to a a given source. 

The source is assumed to be a delta function in space (modelled by a 
sharp Gaussian) with a pulse-like time dependence. The source location, Zg, 
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is on the jet centerline approximately 1*2 jet diameters downstream of the 
nozzle exit# Specifically 


where the function f(t) 


for suitable constants a 
of a broad band spectrum. 


m(t,z,r) 

is 

= f(t)6(r +(z-z ) ), 
s 

(3) 

f(t) 

2 -2 
-(at +bt ) 

= e 

(4) 

and b . 

The use of (4) permits 

the investigation 


The fluctuating field described by (1) reduces to the acoustic field for 
large distances. The near field and flow field are dominated by instability 
waves which are generated because the mean flow is linearly unstable. The 
pulse is assumed to dominate the natural sources of jet noise. These natural 
sources are both the linear and nonlinear terms in (1). In real jet these 
terras are determined from the turbulent fluctuations whereas in the numerical 
model these natural sources are excited by the pulse. The important physical 
effect is the generation of packets of instability waves in the flow. Large 
scale structures which are related to mean flow instabilities have been ob- 
served experimentally in both forced and unforced jets [1,2]. These struc- 
tures interact with and modify the resulting acoustic fields. The model per- 
mits this interaction to be studied in both the linear (e = 0) and nonlinear 
(e ^ 0) regimes. Calculations with the linear model and comparison with 
experiments are described elsewhere [3,4]. In the rest of this paper we 
describe the numerical requirements in order to compute with this model and 
the effects of the nonlinear terms on the fluctuating field. 


Numerical Model 

In order to numerically integrate (1) it is necessary to resolve the so- 
lution over large length scales (far field, near field, and flow field). This 
necessitates the use of a higher order difference scheme. The system can be 
written in the form 


w^ + F + G = H, 
t z r ’ 


(5) 


where w is the vector (p,e,f) and F, G, and H are appropriate 

functions. Equation (5) is split into two one-dimensional operators in z 

and r. Each 1-d system is integrated by using a fourth order version of the 
MacCormack scheme [5] 


w 

i 

n+1 

w. 


“S' * Ik + ''"l 


- >/2 + Is 


At 


( 6 ) 


together with a symmetric variant (Hj^ is obtained from some splitting of 
H) . Typical grids require of the order of 40,000 grid points over distances 
of the order of 50 jet diameters in all directions. Our experience has been 
that on problems of this size, second order schemes are not sufficient to ob- 
tain accurate solutions with the amount of resolution that is feasible. The 
explicit scheme (6) naturally lends Itself to vectorizatlon and has been 
Implemented on the CDC CYBER-203 with great efficiencies. 
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In addition it is necessary to Impose boundary conditions which accurate- 
ly simulate outgoing radiation at the far field boundaries. A family of radia- 
tion boundary conditions has been developed which provide increasingly 

accurate approximations to outgoing radiation. The leading member of this 
family is 


Pt d “ (7) 

w^ere ambient sound speed and the ambient density. Here 

d - r + z and u is the outgoing radial velocity based on a spherical co- 
ordinate system near the source ra* 

System (1) includes terms proportional to r“^. This singularity at the 
axis is resolved by including these terms in the flux vector G when r=0. 
In addition it is necessary to modify the difference formula (6) at 
boundaries. This is accomplished by introducing fictitious grid points 
outside of the computational domain and using a third order extrapolation of 
the flux function (F or G) . This approach was found to be the most readily 
vectorizable. It has been verified that the resulting scheme is fourth order 
accurate* 


Nonlinear Results 

We next describe results illustrating the effect of the nonlinear terms 
on the fluctuating field. In Figs. 2 and 3 the fluctuating pressure is shown 
as a function of axial location z/D (D is the jet diameter) and non-dimen- 
sional time tc^/D for two different radial positions and 
for e - 0.00 and c = 0.05. All figures show an acoustic wave (speed of 
sound normalized to unity) in the downstream direction trailed by several much 
larger waves. These are instability waves which travel with a speed of 
approximately .7Uj where U = .66c^ is the jet exit velocity. A series of 
acoustic ripples can also be seen propagating upstream. These are due to dif- 
fraction of the upstream acoustic wave by the nozzle lip. 

The figures Indicate that the nonlinear terms have little effect on the 
primary acoustic pulse and on the acoustic diffraction from the nozzle lip. 
The nonlinearity has a pronounced effect on the instability waves. Increasing 
the nonlinearity causes these predominantly large scale structures to break up 
into smaller scale structures. This can be seen in both the additional 
ripples which trail the instability waves and a sharpening of the individual 
pulses indicating an enhanced high frequency content. It is also evident that 
for increasing r/D these smaller scale structures are comparable in ampli- 
tude to the primary Instability waves. 

In Figs. 4 and 5 the fluctuating vortlclty is shown for two fixed times 
and for e = 0.00 and e = 0.05. The intense vortices in Fig. 4 correspond 
to the instability waves in Figs. 2 and 3 while the vortices at the later time 
in Fig. 5 represent a residual shedding of vortlclty from the nozzle lip. It 
is apparent from the figures that nonlinear effects tend to slow down the vor- 
tices as they propagate downstream. Thus the trailing vortices catch up with 
the leading vortices and a possible pairing of vortices can be observed. True 
vortex merging, which has been observed experimentally [6], depends heavily on 
viscosity as well as nonlinearity and is not simulated here. 
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In Figs. 6a and 6b the normalized power spectral densities (PSD) for the 
fluctuating axial velocity and pressure are plotted for the linear and non- 
linear cases as a function of the Helmholtz number fD/c^ (f is the 
frequency) . These figures clearly indicate the shift Into higher frequencies 
and the overall broadening of the spectra caused by the nonlinear effects. 
Then results are typical for the fluctuating field at all points. 


Discussion 

The fundamental difference between the nonlinear and linear computations 
is the breakdown of the large scale structures Into smaller scale 
structures. This is associated with a transfer of energy Into higher frequen- 
cies or equivalently a broadening of the spectral content of the fluctuating 
field. It is also illustrated by the increased Interaction between the dif- 
ferent vortices. In real jets this is a fundamental step in the transition to 
fully developed turbulence. The results indicate that at least the Initial 
stages of this energy cascade into smaller scales can be simulated just by the 
nonlinear terms of an axisymmetric and inviscld calculation. 

The generation of smaller scale fluctuations fully justifies the use of 
fourth order spatial discretizations. The numerical scheme is accurate and is 
in general stable. Higher values of e can be readily computed although the 
Gaussian approximation to the 6 function source will have to be smoothed 
out. 
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Figure 1. Computational Domain 
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Figure 3, Three-dimensional plots of the fluctuating pressure r/D = .61 
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5. Fluctuating vorticity at tc„,/D = 30. 





Figure 6a. Flow field velocity spectra z/D = 7.3, r/D = 1.5. 
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